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Abstract 

We measure the masses, decay widths and relative production rate of orbitally 
excited B mesons using 1.25 million hadronic Z decays recorded by the L3 detector. 
B-meson candidates are inclusively reconstructed and combined with charged pions 
produced at the primary event vertex. An excess of events above the expected 
background in the Bit mass spectrum in the region 5.6 — 5.8 GeV is interpreted 
as resulting from the decay B** d — > B^ir, where B** d denotes a mixture of I = 1 
B-meson states containing a u or a d quark. A fit to the mass spectrum yields the 
masses and decay widths of the B^ and Bg spin states, as well as the branching 
fraction for the combination of I = 1 states. In addition, evidence is presented 
for the existence of an excited B-meson state or mixture of states in the region 
5.9 - 6.0 GeV. 
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1 Introduction 



The spectroscopy of orbitally excited B mesons provides important information regarding the 
underlying QCD potential. A flavor-spin symmetry pP arises from the fact that the mass of the 
b quark is large relative to Aqcd- In this approximation, the spin (s*q) of the heavy quark (Q) 
is conserved in production and decay processes independently of the total angular momentum 
(jq = s q + I) of the light quark (q). Excitation energy levels are thus degenerate doublets in 
total spin and can be expressed in terms of the spin-parity of the meson, J p , and the total spin 
of the light quark, j q . The I = mesons, for which j q = 1/2, have two possible spin states: a 
pseudo-scalar P, corresponding to J p = CT, and a vector V, corresponding to J p = l~. If the 
spin of the heavy quark is independently conserved, the relative production rate of these states 
would be V/(V + P) = 0.750 Recent measurements of the B* production rate using B* — > B7 
decays in e + e~ collisions 0IH| agree well with this number. 

For orbitally excited mesons with 1 = 1, there are two sets of degenerate doublets: one 
corresponding to j q = 1/2 and the other to j q = 3/2. For the j q = 1/2 doublet, there is one 
state corresponding to J p = + and three degenerate states corresponding to J p = 1 + . For 
the j q = 3/2 doublet, there are three degenerate states corresponding to J p = 1 + and five 
others corresponding to J p = 2 + . Rules for the decay of the I = 1 states to the I = states are 
dictated by spin-parity conservation [3] . For the dominant two-body decays to a pion and a B 
meson, the j q = 1/2 states undergo an L = transition (S-wave) and their decay widths are 
expected to be broad in comparison to those of the j q = 3/2 states which undergo an L = 2 
transition (D-wave). Table Q presents the nomenclature of the / = 1 B mesons containing 
either a u or a d quark, B** d , with the corresponding spin states, degeneracies and two-body 
decay modes. The spectroscopy of I = 1 B mesons containing an s quark is not studied in this 
analysis, but is examined as a possible source of systematic uncertainty. 





J p 


2J+ 1 


Decay Mode 


Transition 


1/2 


0+ 


1 


B* -> Btt 


S-wave 


1/2 


1+ 


3 


B* -> B*7T 


S-wave 


3/2 


1+ 


3 


Bi -> B*7T 


D-wave 


3/2 


2+ 


5 


B* -> B*tt, Btt 


D-wave 



Table 1: Spin states of the I = 1 B U ( j mesons with the expected relative production 
rates according to spin counting (2 J + 1) and the associated decay modes and 
transitions predicted by spin-parity conservation. 

Predictions for the masses and decay widths of the four spin states are based on Heavy Quark 
Effective Theory (HQET), in which corrections to the spectator quark model are expressed as 
perturbations in powers of Aqcd / "^q [SHE] ■ Such corrections, which can be relatively large for 
c hadrons, are considerably smaller for hadrons containing the more massive b quark. 

Recent analyses at LEP, in which a charged pion produced at the primary event vertex is 
combined with an inclusively reconstructed B meson JH] , have measured an average B** d mass 
in the range 5700 — 5730 MeV. An analysis combining a primary charged pion with a fully 
reconstructed B meson jHj measures Mb* = (5739 (stat) +4 (syst)) MeV by performing a 
fit to the mass spectrum which fixes the mass differences, widths and relative rates of all spin 
states according to the predictions of Reference |H1 

Corrections due to the decay of higher excited states are predicted to be small. 
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The analysis presented here is based on combining a primary charged pion with an inclu- 
sively reconstructed B meson. We use new techniques both to improve the resolution of the 
reconstructed Bit mass spectrum and to unfold this resolution from the signal components. As 
a result, we are able to extract measurements for the masses and widths of both the S-wave 
decays and the D-wave B% decays. 

2 Event Selection and Reconstruction 

2.1 Selection of Z — > bb decays 

The data, collected by the L3 detector [T3] in 1994 and 1995, correspond to an integrated 
luminosity of 90 pb _1 with centre-of-mass energies around the Z mass. Hadronic Z decays 
are selected by making use of their characteristic energy distributions and high multiplicity 
[TT)] . In addition, all events are required to have an event thrust axis direction satisfying 
| cos#| < 0.74, where 9 is the polar angle; to contain a primary vertex reconstructed in three 
dimensions; to contain at least two calorimetric jets, each with energy greater than 10 GeV; 
and to pass stringent detector quality criteria for the vertexing, tracking and calorimetry. A 
total of 1 248 350 events pass this selection. Cutting on a Z — > bb event discriminant based 
on the lifetime information of charged constituents 17J yields a b-enriched sample of 176 980 
events. 

A sample of 6 million simulated hadronic Z decays have been generated with JETSET 7.4 
[T%] and passed through the L3 simulation program [T5] to study the content of the selected 
data. From this study, the purity of the Z — > bb candidates is determined to be 83% and 
corresponds to a selection efficiency of 65%. 

2.2 Selection and Reconstruction 

of b** — > B^tt decays 

The primary event vertex and secondary decay vertices are reconstructed in three dimensions 
on an event-by-event basis such that each charged track can be a constituent of no more than 
one vertex. A calorimetric jet is selected for analysis as a B candidate if it is one of the two most 
energetic jets in the event, if a secondary decay vertex is reconstructed from tracks associated 
with that jet, and if the distance of that vertex with respect to the primary event vertex is 
greater than three times the estimated error of the measurement. 

The decay B** d — > B^7T is a strong interaction and thus occurs at the primary event vertex. 
In addition, the predicted masses for the I = 1 states correspond to relatively small Q values, 
so that the decay pion direction tends to be forward with respect to the B-meson direction 
(Figure P). We take advantage of these decay kinematics by requiring that, for each B-meson 
candidate, there be at least one track which originates from the primary event vertex and which 
is located in the hemisphere defined by the jet thrust axis direction. No attempt is made to 
identify the track as a pion. A total of 60 205 track-jet pairs satisfy these criteria. 

To further decrease background, which is typically due to charged particles from fragmen- 
tation, only the track with the largest component of momentum along the direction of the jet 
is selected. This method has been found [T4*l l2*U] to be an efficient means to improve the purity 
of the signal. In addition, background due to charged pions from D* — > D-zr decays is reduced 
by requiring the track to have a transverse momentum with respect to the jet axis larger than 
100 MeV. These selection criteria are satisfied by 48 022 track-jet pairs with a b-hadron purity 
of 94.2%. 
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The direction of the B candidate is estimated by taking an error-weighted average in the 9 
(polar) and <fi (azimuthal) coordinates of the directions defined by the vertices and by particles 
with a high rapidity relative to the jet axis. A numerical error-propagation method [21] makes 
it possible to obtain accurate estimates for the uncertainty of the angular coordinates measured 
from vertex pairs. These errors, as well as the error for the decay length measurement used in 
the secondary vertex selection, are calculated for each pair of vertices from the associated error 
matrices and determine the weight for the vertex-defined coordinate measurements. 

A second estimate for the direction of the B meson is obtained by summing the momenta 
of all charged and neutral particles (excluding the decay pion candidate) with rapidity y > 1.6 
relative to the jet axis. Such particles have a high probability to be decay products of the 
BW meson. Estimates for the uncertainty of these coordinates are determined from simulated 
B-meson decays as an average value for all events and determine the weight for the rapidity- 
defined coordinate measurements. The final B-meson direction coordinates are taken as the 
error-weighted averages of the two sets of coordinates. 

The resolution for each coordinate is parametrized by a two-Gaussian fit to the difference 
between the reconstructed and generated values. For 9, the two widths are <j\ = 18 mrad and 
02 = 34 mrad, with 68% of the B mesons in the first Gaussian. For 0, the two widths are 
<j\ = 12 mrad and &i = 34 mrad, with 62% of the B mesons in the first Gaussian. 

The energy of the B-meson candidate is estimated by taking advantage of the known centre- 
of-mass energy at LEP to constrain the measured value. The energy of the B meson from this 
method can be expressed as 

j? -fc'cm ~ M-ecoil + jjg f-,\ 

B = 2E ' (1) 

where E cm is the centre-of-mass energy, M reco ii is the mass of all particles in the event with 
rapidity y < 1.6, including the decay pion candidate (regardless of its rapidity) and M B is the 
known B-meson mass. The resolution is estimated by the difference between reconstructed and 
generated B-meson energy values; it is best described by a bifurcated Gaussian [22] with a 
lower width of 1.9 GeV and an upper width of 2.8 GeV. 



3 Analysis of the B7T Mass Spectrum 

The Bit mass is calculated as 



M Bn = ^Ml + ml + 2E B E n - 2p B p n cos a, (2) 

where M B and m n are, respectively, the known B-meson and charged pion masses, E B is the 
B-meson energy described in the previous section, p n is the measured pion momentum and 
a is the measured angle between the B-meson candidate and the decay-pion candidate. The 
resulting mass spectrum is presented in Figure [2K, along with the Monte Carlo background, 
normalized to the region 6.0 — 6.6 GeV. 

The background B7r mass distribution is estimated from the Monte Carlo sample, excluding 
B** — > B^7T decays, and fit with a six-parameter threshold function given by 

Pi(M Bn — P2) P3 e P4 ^ MB ^~ P2 ^ +P5 ^ MB ^~ p ^ 2+P6 ^ MB7T ~ P2 ^ . (3) 

The shape parameters, p^ through p$, are fixed by the fit to the simulated background, while 
the overall normalization factor, pi, is unconstrained in the fit to the data spectrum. In this 
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manner, the normalization uncertainty is accounted for in the statistical error. The uncertainty 
due to the background shape is accounted for in the systematic error estimate. 

To examine the underlying structure, it is necessary to account for effects due to detec- 
tor resolution. The momentum resolution of the decay-pion candidates, which have typical 
momenta of 1 — 3 GeV, is a few percent, and the angular resolution is better than 2 mrad. 
In this case, the dominant sources of uncertainty for the mass measurement are the B-meson 
angular and energy resolutions. Monte Carlo studies confirm that these two components of the 
mass uncertainty are dominant and approximately equal in magnitude. The following analy- 
sis parametrizes the effects of these components on the measured Bit mass and then uses the 
parametrization to fold the resolution effects into the fitting function. 

The dependence of the Bit mass resolution on the B** mass is studied by simulating B** 
decays at four different values of mass and width. Each Btt mass distribution is fit using a Voigt 
function, which is a Breit-Wigner function convoluted with a Gaussian resolution function. The 
Breit-Wigner width is fixed to its generated value and its mass, normalization and Gaussian 
resolution are extracted from the fit. 

Each extracted Bit mass value is found to agree with its generated mass within the statistical 
error of the fit, and the corresponding efficiencies are found to have no mass dependence. The 
Gaussian resolution (Figure Eb) is plotted for each generated B** mass value, and fit with a 
linear function. This function is used to estimate the detector resolution at the measured Bit 
mass. 

Agreement between data and Monte Carlo for the B-meson energy and angular resolution 
is confirmed by analyzing B* — > B7 decays selected from the same sample of B mesons. The 
photon selection for this test is the same as that described in Reference El A B* meson decays 
electromagnetically and thus has a negligible decay width compared to the detector resolution. 
As in the case of the Btt mass resolution, the B-meson energy and angular resolution are 
the dominant components of the reconstructed B7 mass resolution. Fits to the Mb 7 — Mb 
spectra are performed with the combination of a Gaussian signal and the background function 
described above. For simulated events, the Gaussian mean value is found to be Mb 7 — Mb = 
(46.5 ± 0.6) MeV with a width of (11.1 ± 0.7) MeV for an input generator mass difference of 
46.0 MeV. For data, the Gaussian mean value is found to be Mb 7 — Mb = (45.1 ± 0.6) MeV 
with a width of (10.7 ± 0.6) MeV. Good agreement between the widths of the data and Monte 
Carlo signals confirms that the B-meson energy and angular resolutions are well understood 
and simulated. 

According to spin-parity conservation, one expects mass peaks from five possible B** d decay 
modes: Bq — > Btt, B* — > B*tt, B x — > B*tt, B2 — > B*tt and B2 — > Btt. No attempt is made to 
tag subsequent B* — > B7 decays, as the efficiency for selecting the soft photon is low. As a 
consequence, the effective Btt masses for the three decays to B* mesons are shifted down by 
the Mb* — Mb = 46 MeV mass difference. The distribution is fit with five Voigt functions, with 
the relative production fractions determined by spin counting rules. Gaussian convolution of 
the decay widths is taken from the resolution function at each measured mass value. 

Additional physical constraints, presented in Table El are based on predictions common to 
existing HQET models |6TIT2"] . In general, these models predict the mass differences for the two 
j q = 1/2 states and for the two j q = 3/2 states to be approximately equal and in the range 
5 — 20 MeV. Several of the models [310] place the average mass of the j q = 3/2 states above 
that of the j q = 1/2 states, while others [TT1IT2"] predict the opposite "spin-orbit inversion." We 
constrain M B * - M B * = 12 MeV and M B * - M Bl = 12 MeV, but allow the masses of the B{ 
and B2 to be free to test the two opposing hypotheses. 
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k 


Spin State 


Production 


Mass 


Width 


1/2 
1/2 
3/2 
3/2 


B* 
B* 
B, 
B* 


1/12 
3/12 
3/12 
5/12 


M B * - 12 MeV 
free 

M-c* — 1 2 MeV 
free 


r B * 

free 

1 B 2 

free 



Table 2: Constraints applied to the relative production rates, masses and widths of 
the four B** d spin states in the B% mass spectrum fit. 



Predictions are made for the widths of the j q = 3/2 states by extrapolating from measure- 
ments in the D meson system [23! • They are expected to be approximately equal and about 
20 — 25 MeV. No precise predictions exist for the widths of the j q = 1/2 states as there are 
no corresponding measurements in the D system. In general, however, they are also expected 
to be approximately equal, although broader than those of the j q = 3/2 states. We constrain 
Tb* = Tb* and = Tb*, but allow the widths of the B^ and B^ to be free in the fit. 

Simulated signal and background mass spectra are combined and then fit with the functions 
and constraints described above. Mass values and decay widths for the B?; and B^ resonances 
and the overall normalization are extracted from the fit and found to agree well with the 
generated values. All differences lie within the statistical errors. 

4 Fit Results 

The Bit mass spectrum for data is fit with the functions and constraints described above. A 
total of 2770 ± 394 events are in the signal region, and the masses and widths of the B^ and B2 
mesons are found to be 

M B * = (5682 ±23 (stat)) MeV 

T B * = (73 ±44 (stat)) MeV 
M B * = (5771 ± 7 (stat)) MeV 

T B * = (41 ±43 (stat)) MeV. 

The fit has a \ 2 of 81 for 70 degrees of freedom. It does not describe the data well in the 
region 5.9 — 6.0 GeV, where there appears to be an excess of data events over the simulated 
background. Exclusion of this region from the fit yields consistent values for the signal and 
reduces the x 2 to 57 for 60 degrees of freedom. 

Several HQET models [7[ l9* [IT2^ predict the existence of radially excited (2S) B-meson states 
in the region 5.9 — 6.0 GeV. In addition, an inclusive measurement of Btctt final states in Z 
decays [21] provides evidence for a resonance in the same mass region. To account for the 
possible existence of these states in the mass spectrum, we refit the data including a Gaussian 
function in the region of interest. 

The resulting fit, shown in Figure [31 has a \ 2 °f 63 for 67 degrees of freedom. A total of 
2784 ± 274 events are in the original signal region, corresponding to the branching fraction 
Br(b -> B;* d -> BWtt) = 0.32 ± 0.03. The masses and widths of the BJ and B^ mesons are 
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found to be 




(5670 ± 10 (stat)) MeV 
(70 ± 21 (stat)) MeV 
(5768 ± 5 (stat)) MeV 
(24 ± 19 (stat)) MeV. 



In addition, a total of 297±100 events are in the high-mass Gaussian, denoted B', corresponding 
to the branching fraction Br(b — ► B' — ► Bit) = 0.034 ± 0.011. The mass and Gaussian width 
are found to be M w = (5936 ± 22) MeV and a B > = (50 ± 23) MeV. Table El presents the 
correlations of all free parameters in the fit. As this fit yields the best overall confidence level 
(61.2%), compared to the other two fits described above, its values are chosen as the results of 
the analysis. For comparison in the studies below, we refer to it as the "final fit." 



Param. 




iV B ** 


M B * 




M B * 




N B , 


M w 


cr B ' 


^bgd 


1.000 


-0.579 


0.144 


-0.140 


0.100 


-0.231 


-0.606 


-0.157 


-0.432 




-0.579 


1.000 


0.004 


0.468 


0.133 


0.607 


0.050 


0.289 


-0.022 


M B * 


0.144 


0.004 


1.000 


0.023 


0.276 


0.269 


-0.244 


0.090 


-0.190 




-0.140 


0.468 


0.023 


1.000 


-0.124 


0.506 


-0.115 


0.120 


-0.086 


M B * 


0.100 


0.133 


0.276 


-0.124 


1.000 


-0.041 


-0.355 


0.247 


-0.362 


r B ^ 


-0.231 


0.607 


0.269 


0.506 


-0.041 


1.000 


-0.127 


0.192 


-0.060 


N w 


-0.606 


0.050 


-0.244 


-0.115 


-0.355 


-0.127 


1.000 


-0.024 


0.773 


M w 


-0.157 


0.289 


0.090 


0.120 


0.247 


0.192 


-0.024 


1.000 


0.058 




-0.432 


-0.022 


-0.190 


-0.086 


-0.362 


-0.060 


0.773 


0.058 


1.000 



Table 3: Table of correlations of all free parameters in the final fit. 

The results of the final fit place the average mass of the j q = 3/2 states (98 ±11) MeV 
higher than that of the j q = 1/2 states. This supports some theoretical predictions [3121 but is 
contrary to those which predict spin-orbit inversion [THIT2"] . To test the ability of the procedure 
to discriminate between the two possible cases, additional fits are performed. In these fits, the 
widths and relative production rates of the j q = 1/2 and j q = 3/2 states are varied, while the 
masses of the j q = 1/2 states are constrained to be either equal to those of the j q = 3/2 states 
or up to 100 MeV higher. In all cases, the fit confidence levels decrease as the mass of the 
jq = 1/2 states is increased relative to that of the j q = 3/2 states. The highest confidence level 
for a fit supporting spin-orbit inversion is found to be 8.3%, relative to 61.2% for the final fit. 

The fits performed above constrain both Mb* — M Bl and Mb* — Mb* to 12 MeV. The 
results of varying this value to test the constraint as a source of systematic uncertainty are 
presented below. It was suggested to independently test the validity of the j q = 1/2 mass- 
difference constraint since no corresponding measurements exist for the D-meson system, as 
they do for the j q = 3/2 states. This test is performed by fitting the mass spectrum with 
Mb* — Mb* as an additional free parameter. The resulting fit, which has a x 2 °f 63 for 66 
degrees of freedom, yields masses, widths and branching fractions for the B^, BJ; and B' which 
are consistent with the final fit. The mass difference of the j q = 1/2 states is measured to be 
Mb* — Mb* = (—39 ± 71) MeV. While this value is consistent with the original constraint, the 
large error indicates that the method is insensitive to the mass difference of the j q — 1/2 states 
and, conversely, that the results of the final fit are only weakly dependent on that constraint. 
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As an additional check, the mass spectrum is fit with the combination of a single Voigt 
function and the background function. The resulting fit has a \ 2 °f H4 for 76 degrees of 
freedom and describes the signal region poorly. A total of 1854 ± 153 events are in the signal, 
denoted B**, corresponding to the branching fraction Br(b — > B** — > B^n) = 0.21 ± 0.02. 
The mass and width are found to be M B ** = (5713 ± 2) MeV and r B » = (31 ± 7) MeV. The 
hypothesis that the signal be the result of the decay of a single resonance is highly unlikely 
considering the low confidence level of the fit (0.35%). 

Sources of systematic uncertainty and their estimated contributions to the errors of the 
measured values are summarized in Table EJ The b-hadron purity is varied from 92% to 
96%. This variation affects only the overall B** and B' production branching fractions. The 
contribution to these errors due to uncertainty in R^, the branching ratio Br(Z — > bb)/Br(Z — > 
qq), is negligible. 



Sources of 


SM Bt 


5T B - t 


SM B * 


5T B * 


^Br(B-) 


SM B / 


5a B i 


<SBr(B') 


Uncertainty 


[MeV] 


[MeV] 


[MeV] 


[MeV] 




[MeV] 


[MeV] 




b purity 










±0.01 






±0.007 


background 


±3 


±19 


±2 


±12 


±0.06 


±2 


±2 


±0.002 


mass constraints 


±4 


±7 


±4 


±7 


< 0.01 


±1 


±1 


< 0.001 


width constraints 


< 1 


±2 


< 1 


±2 


< 0.01 


< 1 


±2 


< 0.001 


fraction B* -> B, B* 


±1 


±2 


±2 


±4 


< 0.01 


±3 


±4 


±0.003 


signal resolution 


±1 


±9 


±2 


±9 


±0.01 


< 1 


< 1 


< 0.001 


signal efficiency 


< 1 


±1 


< 1 


±1 


±0.01 


< 1 


< 1 


±0.001 


inclusion of B' 


±12 


±12 


±3 


±17 


±0.01 








Total 


±13 


±25 


±6 


±24 


±0.06 


±4 


±5 


±0.008 



Table 4: Sources of systematic uncertainty and their estimated contributions to the 
errors of the measured values. Effects due to correlations between the quantities are 
taken into consideration in the total error. 



Systematic effects due to background modelling are studied by varying the shape parameters 
of the background function (parameters P2 through p 6 ) and by performing the fit with other 
background functions. Effects due to signal modelling are studied by varying the constraints 
on the masses, widths and relative production rates. The mass differences M B * — M Bl and 
Mb* — Mb* are each varied within the range 5 — 20 MeV, covering the full range of theoretical 
predictions. The ratios of the Breit-Wigner widths I^/rB* and TB*/rB* are each varied 
between 0.8 and 1.0. This is to account for the possibility of a slight mass dependence due to 
the available phase space and the probability for B** — > B^p decays Finally, the fraction 
of Bg mesons decaying to B* mesons is varied between 1/3 and 2/3. 

Effects due to uncertainty in the resolution function are estimated by varying the slope and 
offset of the linear parametrization. Effects due to mass-dependent variations in the signal 
efficiency are estimated by performing a linear parametrization of the Monte Carlo efficiencies, 
similar to that of the resolution, and then varying the corresponding slope and offset. 

Bit pairs from the decay B** — > Birir, for which only one of the pions is tagged, are studied 
as a possible source of resonant background. The resulting reflection is found to contribute 
insignificantly to the background in the low Q-value region. Similarly, B** — > BK decays, for 
which the K is misidentified as a 7T, are found to contribute only slightly to the low Q-value 
region, and their effects are included in the background modelling uncertainty contribution. 
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5 Conclusion 



Inclusively reconstructed B mesons in Z — > bb events are combined with charged pions orig- 
inating from the primary event vertex. A fit is performed to the B7r mass spectrum in the 
framework of Heavy Quark Symmetry pQ and under constraints common to several theoretical 
models. The results of the fit, which provide the first measurements of the masses and decay 
widths of the (j q = 1/2) and B2 (j q = 3/2) mesons, as well as the branching fraction for the 
combination of I = 1 states, are 

(5670 ± 10 (stat) ± 13 (syst)) MeV 
(70 ± 21 (stat) ± 25 (syst)) MeV 
(5768 ± 5 (stat) ± 6 (syst)) MeV 
(24 ± 19 (stat) ± 24 (syst)) MeV 
0.32 ± 0.03 (stat) ± 0.06 (syst). 

In addition, an excess of events near 5.9 — 6.0 GeV is interpreted as resulting from the decay 
B' — > B^*^7r, where B' denotes an excited B-meson state or mixture of states. From the same 
fit, the B' mass, Gaussian width and branching fraction are 

M w = (5937 ± 21 (stat) ± 4 (syst)) MeV 
a B > = (50 ± 22 (stat) ± 5 (syst)) MeV 
Br(b -> B' -> B w 7r) = 0.034 ± 0.011 (stat) ± 0.008 (syst). 

For both branching fractions, isospin symmetry is employed to account for decays to neutral 
pions. 

These results disfavor recent theoretical models proposing spin-orbit inversion JUE], but 
agree well with earlier models [3113 an d provide support for Heavy Quark Effective Theory. 
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r B * = 

Br(b -> B** d -> B w tt) = 
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Figure 1: Schematic of the decay B** — > B^ir. Ellipses represent the reconstruction 
uncertainties of the vertex positions. The decay pion points to the primary vertex 
and its direction is forward with respect to the direction of the B meson. 
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Figure 2: (a) Mass spectrum for selected Btt pairs. The points with error bars 
are data and the histogram represents the expected background from Monte Carlo, 
normalized to the sideband region 6.0 — 6.6 GeV. (b) Linear fit of the extracted Bn 
mass resolution for the Monte Carlo signal components at the generated B** mass 
values. 
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Figure 3: (a) Fit to the data Bit mass distribution with the five-peak B** signal 
function, the Gaussian B' signal function and the background function described in 
the text, (b) The resulting background-subtracted distribution. 
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